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I. INTRODUCTION 



The angular part of the Faddeev equations is solved analytically for s-states for two-body square- 
well potentials. The results are, still analytically, generalized to arbitrary short-range potentials for 
both small and large distances. We consider systems with three identical bosons, three non-identical 
particles and two identical spin- 1/2 fermions plus a third particle with arbitrary spin. The angular 
ly-^ ' wave functions are in general linear combinations of trigonometric and exponential functions. The 

Efimov conditions are obtained at large distances. General properties and applications to arbitrary 
potentials are discussed. Gaussian potentials are used for illustrations. The results are useful for 
numerical calculations, where for example large distances can be treated analytically and matched 
^ I to the numerical solutions at smaller distances. The saving is substantial. 

'nI" ! PACS numbers 3.65.Ge, 2.60. Nm, ll.80.Jy 21.45. -fv 
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new method to study the quantum mechanical three-body problem was recently formulated . The first signifi- 
. cant advantage is the precise treatment of large distances where the correct asymptotic behavior of the wave function 
is incorporated. Another advantage is the parallel treatment of bound and continuum states, which in combination 
with the long-range treatment opens the interesting possibility of approaching the celebrated three-body Coulomb 
^ problem in the continuum. The method solves the coordinate space Faddeev equations in two steps. First the angular 
f5 i (and most difficult) parts of the three-component wave functions are calculated by elaborate use of the analytical 
knowledge of the large-distance behavior. Then the coupled set of effective radial equations is solved numerically [||. 
• ^ The new method is designed to solve the Faddeev equations in coordinate space with the advantage of an immediate 

intuitive interpretation of the physics involved. (It may be considered as a development from the approach used to 
study the properties of [^.) The advantages are seen in the analytical and numerical treatment of extreme cases 
like the so-called Efimov states , which has been suggested and looked for in both molecules ||] and nuclei 1^ . These 
states occur when at least two of the two-body subsystems simultaneously have an s-state at zero energy arising from 
short-range interactions. The resulting infinitely many bound three-body states of 0'*' nature are extremely extended 
in space and extremely weakly bound. The long range Coulomb interaction between only one pair of particles destroys 
the effect. 

The large-distance coupling takes place only between s-states in different relative two-body subsystems ||^ . Further- 
more, low-lying three-body states often contain large components of such s-states. Different treatments are usually 
needed when long-range interactions like the Coulomb potential are involved and we shall only consider short-range 
interactions. We shall furthermore neglect the coupling to higher angular momenta and confine ourselves to relative 
s-states. The purpose of this paper is to classify the lowest eigenvalues and describe how the solutions for small as 
well as large distances can be obtained analytically for the angular part of the Faddeev equations. These solutions 
turn out to be exact for square-well potentials where we also can find exact solutions at intermediate distances. Thus 
we shall derive a semianalytic s-state square-well solution to the Faddeev equations. 

The solutions describe the Efimov anomaly and how the infinitely many states continuously appear and disappear 
as function of the parameters of the potential. The power of the method is especially due to the analytical treatment 
of large distances which must be treated with particular care for loosely bound quantum systems in low angular 
momentum states. Also the connection between two- and three-body large-distance behavior can be explored. Using 
the analytical results for small and large distances in numerical calculations for arbitrary potentials improve both 
precision and computational speed and enable thereby investigations of otherwise inaccessible problems. 
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Other methods are available for investigations of the three-body problem for baryons 0, for Borromean systems 
for molecular systems as well as for Coulomb interacting particles The results from all these 

methods are needed in studies of the properties of a variety of different systems of interest in physics, see for example 
-pTf. Some of the methods solve the Schrodinger equation directly, but the Faddeev equations are needed to 
describe all the subtle correlations. This is especially clearly seen in systems where two- and three-body asymptotic 
behavior are mixed Borromean systems, characterized by genuine three-body asymptotics only, are relatively easy 
to handle, whereas weakly bound systems, where both two and three-particle correlations are essential, require very 
careful treatment of the large distances. 

Analytic results are rare in quantum mechanics where the Schrodinger equation should be solved. The Faddeev 
equations further complicate analytic analyses. However, so far at least one exception exists for identical spinless 
particles interacting via two-body harmonic oscillator potentials . This potential is infinite at large distance and 
scattering states therefore cannot be studied. Also the behavior for short-range potentials in general are excluded. 

In this paper we outline in section 2 the general theoretical framework and in section 3 we solve the angular 
eigenvalue problem for s-states for a system of three identical bosons and spin-independent interactions. In section 
4 we generalize to systems of three different particles and in section 5 we consider a system with two identical spin- 
1/2 particles plus a third particle. In section 6 we give numerical illustrations and indicate qualitatively how to 
generalize the results. Tedious mathematical derivations are collected in appendices. Finally we give a summary and 
the conclusions in section 7. 



II. THEORETICAL FRAMEWORK 



The intrinsic Hamiltonian of the three-body system is given by 

i—l i>j — l 

where m.^, r.^ and are mass, coordinate and momentum of the z'th particle, Vij are the two-body potentials, P and M 
are respectively the total momentum and the total mass of the system. We shall use the (three sets of) hyperspherical 
coordinates which consist of one radial coordinate p (hyperradius) and five generalized angles Qi, where i= 1, 2, 3. The 
precise definitions are given in appendix A. One of these sets of hyperspherical coordinates is in principle sufficient 
for a complete description. The volume element is given by p^dfldp where dfl — sin^ acos^ adadQ^d^y 



A. General procedure 

The total wave function of the three-body system is written as a sum of three components ■0^*'' which in turn for 
each p are expanded in a complete set of generalized angular functions: 

* = EV'^'^ = , (2) 

2—1 ^ i,n 

where the radial expansion coefficients fn{p) are component independent and p^'°l'^ is the phase-space factor. These 
wave functions satisfy the three Faddeev equations 

(T - -f VSfe(^^*^ -f + ^^''^) = , (3) 

where E is the total energy, T is the kinetic energy operator and {«,_;/, /c} is a cyclic permutation of {1,2,3}. The 
Faddeev equations may have non-trivial spurious solutions, where each component is non-vanishing while the sum 
corresponding to the Schrodinger wave function is identically equal to zero. The components of such solutions are 
eigenfunctions of the kinetic energy operator with eigenvalues equal to the total energy. 

The procedure is now for each p to solve the eigenvalue problem for the five dimensional angular part of the Faddeev 
operator: 

|1 1 A^<i>W + y^^(^ + $0) + ^ £ 1 A„(p)<i>W , (4) 
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where is the p-independent part of the kinetic energy operator defined by 



Exphcitly the generahzed angular momentum operator is given by 

1 P P 

A^ = — : — ^ sin a cos a H % 1 \ 4 (6) 

sm a cos a aa^ sin a cos^ a 

in terms of an arbitrary a-coordinate and the angular momentum operators l"^ and related to the Jacobi coordinates. 
The spurious states are characterized by angular eigenvalues equal to those of the angular kinetic energy operator A^, 
i.e. K[K + A),K = 0,1,2... 

Insertion of tp^"^^ defined in cq. (||) into eq. (||) then leads to the coupled set of "radial" differential equations 

^ - ^ + ^ (a„(p) + f )) = E {^Pnn'^^ + Q..) , (7) 

with the functions P and Q defined by 



i,3 = i 



(8) 




The diagonal part of the P-matrix vanishes, i.e. Pnn = 0. 



B. Angular eigenvalue equation 



It is convenient to show explicitly the spin dependence of the wave function ^n\p, ^i) in cq. (^, see e.g. Let 



us consider s-waves only and assume that s is the intermediate spin obtained by coupling of the spins Sj and of 
particles j and k. The spin s is afterwards coupled to to give the total spin S of the system. The total angular 

t>n^s{p,^i) and the spin dependent part Xn,s, 



wave function for the i'th channel then factorizes into the spatial part 
i.e. 



^iHP^^^ 



1 



sin {2ai 



(10) 



where we explicitly extracted the phase-space factor sin (2cki). Both (j) and x may depend on the intermediate coupling. 
The Faddeev components in eq. (|^) must be expressed in one Jacobi coordinate set. For the s-waves the wave 

(k) 

functions (j)n,sj which only depend on ak and p, are first expressed in terms of the i'th set of hyperspherical coordinates 
and subsequently integrated over the four angular variables describing the directions of and y^. This amounts, 
still for s-waves, to the substitution formally expressed by the operator Rik defined by 



Ak) 



sin(2afe) 



shi{2(pj) sin(2ai) J\^. 



TT/2-\n/2-ipj-ai\ 



4^liP,ak)dak , 



(11) 



where the angle ipj is given by the masses as 



tan ipj 



I rrij {mi + ni2 + ttis) 



mitrLk 

It is closely related to the transformation angle ipik used in appendix A and defined by 



(12) 
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'fik — — arctan 



mj{mi +1712 + mz) 
y mivrik 



(13) 



where (—1)^ is the parity of the permutation {i, fc, j}. 

Substituting eq. into eq. (Q) we obtain, after multiplication from the left with Xn]s, the angular eigenvalue 
equation 



da? 



+ {p''{xl:\W.{psma^)\xl:b - \n{p)WAP, a^) + p2sin(2a,) 







r aC') 1 












sin(2aj) 


sin(2Qfe) 



(14) 



where Vi{x) = Vjk{x/ajk)'2m/h with ajk defined in appendix A, A„(/9) ~ A„(/9) + 4 and the coefhcients C*^, expressing 
the overlap of the spin functions are given by 



Qlk ^ / (i) I W \ 



(15) 



These matrix elements are diagonal for i — k, i.e. C*^, = Sgs' and symmetric, i.e. Cf^, = C^^. 
When the potential is diagonal in spin-space we obtain the much simpler set of equations 



da? 



+ (p'(xEl^'.(psinaOlxE> " A„(p))0l*Up, + p'sin(2a0 



:(x«|«.(psinaO|x«)EK^-'^^. 



sin(2aj) 



r^ik r? 
^ss'^ik 



sin(2afe) 



(16) 



Eqs. ( PI ) and (|6|) constitute sets of equations obtained for i=l,2,3 and all possible (i-dependent) values of s. As 
usual the values of {«, j, k} must here be a permutation of {1, 2, 3}. 



C. Spin independent interactions 

When the interactions are independent of spin, each Faddeev component must factorize into a spin-part and a 
spatial part. Furthermore, the spin can be factorizcd out of eq. (^) and the structure of the Faddeev equations then 
remains unchanged for the spatial parts alone. The spin dependent wave function must then be the same for all three 
Faddeev components and the spatial parts of the wave function in eq. ( p^ ) can at most differ by a normalization 
constant, i.e. 



0«(p,l])^&W0«(p,f]), 
E^.^.^E^i'^X^'^ z,fc=l,2,3. 

S S 

Multiplication of eq. ( p^ ) by Xn,s from the left then gives 

s' 

These equations are not independent, since 

s' s' 

for all values of 5, k and I. This is easily seen by use of the closure relation of eq. (fq), i.e. 



(17) 
(18) 

(19) 
(20) 
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ss' — 2^^ss"^s"s' •> t — 1,Z,0 



(21) 



Thus eq. (|l^) is valid for all k, if it only holds for one of the values of k. If b'^f' is arbitrarily chosen and b\''' and 
foi^'' calculated (i=l,2 and k=3) from eq. (^9|), we can see by using eqs. (|2l| ) and (^o|), that eq. (^9|) is valid for all 
other values of i and k. Thus any choice of llt^ for one value of i provides the same spin-independent solution of the 
Faddeev equations. 

In the symmetric case when the three particles furthermore have equal masses and spatial interactions, the wave 
functions $1*'' corresponding to the different Faddeev components are independent of i. Then eq. (10) implies that 



^ (p, 0)X« = ^n}s (P, ^)xi'l , 



(22) 



which by use of eq. ( |17[ ) directly leads to eq. ( |18D 

Inserting eq. (|l^) into eq. ([l^) we obtain by use of eq. (|l^) the equation 



d'^(t>n{p, aj) 



which determines 



+ p sin(2Q!i)ui(psinai) yRij 
= ch'i^ for i=l,2,3. 



sin(2aj) 



Rik 



_sin(2afc) 



0, 



(23) 



III. THREE IDENTICAL BOSONS FOR SPIN INDEPENDENT INTERACTIONS 

In the symmetric case with spin independent interactions the remaining single Faddeev equation eq. (|2^ ) can by 
use of eqs. ( |ll| ) and (^2|) explicitly be written as 



d'^(f>ip,a) 



+ (p^w(psina) — X{p))4>{p, a) 



—/=p w(psina) 
v3 



ir/2-|ir/6-Q| 



|ir/3-Q| 



(/)(p, a')da' — , 



(24) 



where we omitted the label "n" . We shall only consider short-range potentials and often further restrict ourselves 
to square wells. It is most convenient first to solve exactly for the schematic square-well potential and afterwards 
generalize as much as possible. The two-body potential is then a step function 



and the reduced potential in eq. 



V{r) = -VoSir < Ro) 
is therefore another stcpfunction 
v{x) = -vqQ{x < X) , 



(25) 



(26) 



where vq = 2mVo/h^ and X — Ro/\/2 when the normalization mass is equal to the mass of one the particles. The 
angular Faddeev equation in eq. ( p4|) is now solved analytically in different intervals corresponding to increasing values 
of p. The decisive quantity for the square- well potential is psina, which determines whether v is finite or vanishes. 
When p < X we have v = —vq for all values of a e [0, n/2]. 



A. Short-distance behavior: < p < Ro/V^ 
The potential is now constant for all a and eq. (pj) simplifies to 



<9 0(p,a) .2 , w N 4 2 

-{p vo + \{p))cj){p, a) = — =p uo 



/2-|7r/6-a| 



(/)(p, a')da' 



(27) 
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The differential equation, without tlic right hand side, has vanishing solutions at a = of the form sin(fca), which by 
insertion on the right hand side of eq. ( p7|) e asily is shown to remain of this form if k is an even integer, i.e. k — 2n, 
where n is an integer. A solution to eq. (|27D, 

(f>{p, a) cx sin(2nQ;) , (28) 
is then by insertion found to correspond to the values of An given by 

A„ = 4n2-p2„o(^l-^(-)"^sin(n7r/3)^ . (29) 

This wave function is furthermore also the solution for all potentials at p = where the corresponding eigenvalue 
Xn = 4,n'^ then is obtained by demanding a vanishing wave function at a = 7r/2. Since the potential only enters in 
combination with p^, a perturbative solution to second order in p^ is easily obtained for potentials, which are finite 
at the origin, by using v{x) « w(0) ~ — vo- Thus eqs. ( |29| ) and (^) are also solutions for any potential in first order 
perturbation theory provided the depth of the square- well potential is replaced by the value v(0) of the potential at 
P^O. ' 

The solution for n = 2 is independent of p and corresponds to the trivial solution (identically zero) to the Schrodinger 
equation. For n > 2 the solution is only an approximation, since the higher angular momenta (neglected here) 
contribute to these states. 



B. Intermediate distances: 7r/3 < ao < i" /2,Ro/V2<p< Roy/2/3 
The potential now vanishes for ao < a < tt/2 (region II) where 

ao = arcsin(i?o//OV^) (30) 

and it remains constant and finite for < a < ag. The solution to eq. (|2j) in region II (vanishing potential) is 
therefore 

077 = An Sin ((a - 7:/2)Wj = ^^77 (^e*("--/2) VI _ ^~^{a^./2)Vl^ ^ (3^) 

where we already explicitly selected the solution vanishing at a = 7r/2. 
To proceed we divide the a-space into subregions: 

A: 0<a<ao-^ 

B:ao~-<a<—-ao (32) 

^ 27r 

: - "0 < a < "0 

// : ao < a < 7r/2 . 

These regions are marked in fig. 1, which shows the integration limits of a' as function of a for the integral in 
eq. (p4[). The coupling between the different regions can then be seen and expressed in the set of equations: 

- (p vo + X(p))(pA(P,a) ^ —^p vq (l)c{p,a)da , (33) 



|7r/3-a| 



^^^^^ ~ {p'v, + \{p))^c{p. a) = ^yv, (^7/3:/; 0^(p, a')da' 

+ &3Z ^c{p, a')da' + C:J-;,T Mp, «')d«') , (34) 
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71/2 
71/3 

a' 

27u/3-a„ 



71/6 



a.g-Ti/3 












ao-Tz/3 7x;/6 27i;/3-ao 

a 



71/3 



71/2 



FIG. 1. The integration area in eq. ( p4[ ) is the tilted rectangle (sohd hnes). The regions A,B,C and II (see eq. (p2[)) are 
marked on the figure. The size of ao corresponds to intermediate values of p: Ro/V^ ^ P ^ i?o-\/2/3. The thick solid line on 
the axis indicates the range of ao— values. 



doL^ 



{p^vo + \{p))4>B{p-,a) = ■^p'^vq (^J^J/^_"\° 'i>B{p,a')da' 



(35) 



where the solutions are labeled according to subregion. Only (p^ and (pc are directly coupled since 4>b enters in 
eq. ( p4[ ) as an integral over a constant interval. Analogously only pB and pu are directly coupled. 

Integrating the (^//-term in eq. ( ^5| ) we obtain exponentially increasing and decreasing functions, exp(±ia\/X), 
which by further integration and differential derivation still remain of the same functional form. Thus such functions 
matching (/>// are necessary in the solution, but in addition other exponentials are also possible. By insertion we then 
find that the wave function 



IIia\/X 



fc=i 



B^_^K,- 



(36) 



is a solution to eq. (pa) when the B-coefficients and 



are related as shown in appendix B. We have unique 
solutions for B^^ and B[/ , except the patological case of A = 16/3, and three different solutions for k^-* with three 



corresponding constraints between s'j*'' and B^' . This is the explanation for the three terms of the same form in 

eq. (^6|). In addition there is a link to the C-region providing one constraint between pc and the _Bj_'^ ''-values. 

In the coupled A- and C-regions we must also look for exponential functions as solutions. By insertion into eqs. (|32 
and (^) we find then that the wave functions 



k=0 



AC 



(37) 



3 

E 

fc=0 



(C. 



{k)^aK 



(k) 

AC 



(38) 
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indeed are solutions provided the A- and C-coefRcients and the K^^-values are related as shown in appendix C. 
We have in eq. (5^) already imposed the constraint that (j)A{a = 0) = which eliminated one of the coefficients 

(k) 

present in the other regions. We obtain four different solutions for as indicated by the four k- values and the two 
C^'^^-coefficients are uniquely determined by A^''\ In addition there is a link to the B-region providing one constraint 
between (pB and the A'^'^^-values. 

The solutions are now explicitly written down in the four regions named A,B,C and II. We have found wave functions 
containing the 8 parameters An, b'^\ A^^^ and A^'^'' for k=l,2,3. Two linear constraints exist between them as seen 
from eqs. (b|1i|) and (b|T^) in appendix B and eqs. (C||) and (C|lO|) in appendix C. The matching conditions at the 
three boundaries between the regions then provide additional 6 linear constraints on the remaining 6 free parameters. 
This leads as usual to the quantization condition for the eigenvalue A. 



n/2 



a 



n/3 



71/6 



7r/3-a„ 








7C/2 



FIG. 2. The same as figure 1 for an ao— value corresponding to intermediate values of p: _Ro\/2/3 < p < Roy/2. The regions 
marked on the figure are defined in (see eq. (P9[)). The thick solid line on the axis indicates the range of ao— values. 



C. Intermediate distances: 7r/6 < ao < tt /3 , Ro Jiji <p<RoV2 



The potential again vanishes for ao < a < tt/2, where ao is given in eq. (|30|). Subregion C has now been absorped 
into region II and subregions A and B are, as shown in fig. 2, now defined by 



The coupled set of equations are now 



A: 0<a< ao 

- - 3 

TT 

B ■■ - - ao < a < ao 
J/ : ao < a < 7r/2 . 



(39) 



4 

(p'^vq + X{p))(j)Aip, a) = —r^P^vo I (l^iiiP, a')da' , 

V3 Jtt/S-q 



(40) 



d <pB{p,a) 



{p^vo + \{p))(t)B{p, a) ^ -^P^vo (/"/Va ^b{p, a')da' 
+ Jao (l)ii{p,a')da' \ . 



(41) 
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The couplings are much simpler and essentially only (pg enters in an integro-differential equation. Since 0// still is 
given by the expression in eq. (^l|), a solution to cq. ( |40| ) would have to be proportional to sui[a\/\) . The solution, 
where the right hand side vanishes, is analogously proportional to sin(Q!K), where 



K = \l p^vq + \{p) 

By insertion we then obtain the complete solution to eq. (Eo[) as 



'A _ g-«"> 



1 



4) A = Af{e 

where the coefficient Ah is arbitrary and ^/ is given by 

').Att . P, 



+ -Ahie'""' - e^'""") = 2iAf sin 



m + Ah sin(Q;K) 

in(^aVA7r/6) . 



(42) 
(43) 
(44) 



/3A V3A 
The wave function in the B-region is a solution to eq. ( p|) , which apart from a constant term arising from if)c and 
the upper limit of the integral of ipB, is identical to eq. (35). The solution is therefore given by eq. (pq ), which by 
insertion into eq. ( |4l| ) also in this case leads to the expressions for the coefficients given in eqs. (B|i|) — (B^in appendix 
B. The only difference is the constants in eq. (B|), which now is changed into 



kl. ( gi\/I(Qo-7r/2) ^ g-i\/J(Qo-7r/2)\ ^ 1_ / 

+ J2k=i ( 



B 



(45) 



With the expressions in eqs. (B§), (B0) and for the coefficients we can rewrite the constraint from eq. (|4^) as 



(7r/3-ao) _ g-K^*'(27r/3-Qo) 



+ 



2VA 



^gjAA(7r/2-Qo) _|_ g-i-\/I'(7r/2-Qo)) 



(46) 



giVA(Qo-7r/6) 2 - 4 ^ 



g -jVA(ao-^/6) [ 2 + 4 ^ 



a' 




«o 71/6 71/3 71/2 

a 

FIG. 3. The same as figure 1 for an ao— value corresponding to large values of p: RoV^ < p < oo. The regions marked on 
the figure are defined in (see eq. (^)). The thick solid line on the axis indicates the range of qq— values. 
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We are now left with 5 parameters, i.e. Aji, Ah and si*^-* for k=l,2,3. One of these is determined by the constraint 
in eq. ( ^ ) . The matching conditions at the two boundaries between the regions provide additional 4 linear constraints 
on the remaining 4 free parameters. Therefore we obtain again a quantization condition for the eigenvalue A. 



D. Large-distance behavior: < qq < 7r/6 , RoV^ l£ P l£ oo 



The potential vanishes as usual for ao < a < tt/2, where ao again is given in eq. (pOj). Subregion B has now also 
been absorped into region II and subregion A is, as shown in fig. 3, now defined by 



A: < a < an 
II :ao<a< 7r/2 



(47) 



The equations now reduce to eq. (pO|), which couples (pA and </>//. The latter is again given by the expression in 
eq. ( 31 ) and consequently (pA is given by eq. ( ^3|) with one arbitrary coefficient and the other coefficient expressed in 
eq. (Q. 

The matching conditions at ao provide 2 linear constraints on the remaining 2 free parameter and we obtain again 
a quantization condition for the eigenvalue A. An explicit expression is most conveniently derived by rewriting the 
solution in eq. (El) as 



8A 



II 



= sin(v A 7r/6) sin(aV A) + Ah sin(Q!K) 



Matching the wave functions in eqs. ( |3T| ) and (|4^) and their first derivatives at a = qq gives immediately 



8A 



'3A 



= sin('\/X 7r/6) sin(ao'\/X) + Ah sin(aoK) = Ajj sin ((ao — 7r/2)'\/X) 



(48) 



(49) 



V3 



sm 



A 7r/6) cos(aoVA) + AhK,cos{aQK) = Ajj\/\cos (j^a^ — tt/2)\/\ 



The determinant D corresponding to this linear set of equations for An and Ah is 



D = -j^ cos(aoK) 



V3' 



A 7r/6) sin(ao\/A) — ■\/Asin ^(ao — tt/2) 



— sin(aoK) sin(\/A n/6) cos(ao\/A) — VTcos (^{ao — n/2)\/\^ 

The non-trivial solutions are obtained for D = 0, which therefore is the eigenvalue equation for A. 
We can simplify even further in the limit of very large distances, where 



ao 



Rq 



(50) 



(51) 



(52) 



where we assumed that A remains finite when p increases towards infinity. The eigenvalue equation D = obtained 
from eq. (|5l|) by expansion to the lowest orders in 1/p then becomes 



—= sin(\/X 7r/6) — \/\cos{\/\ tt/2) 

v3 O'scat 



sm( V A tt/2) 



where we introduced the two-body scattering length ascat of the square-well potential 

1 VV2 



ascat -~Roy/vo/2 + tan(i?o\/wo/2) 



(53) 



(54) 



The factor of two appearing here is due to the mass (not the reduced mass) in the definition of vo ■ When p approaches 
infinity eq. ( |53| ) can be valid only if A — > An^, where n must be an integer. This is the same spectrum as obtained in 
eq. ( p9|) for p ^ 0. Expanding to the next order in 1/p gives 
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( 1 



T^pV2 ) 



(55) 



If the scattering length is infinitly large eq. ( p3| ) reduces to the equation derived by Efimov Q| , which instead has 
the solution A = —1.01251 leading to infinitly many bound states obtained from eq. (^. Clearly if need be these 
expressions can be extended to any order of correction in 1/ p. 

When X/ diverges the asymptotic behavior is different. We then define in the large-distance limit by A 
2mtxp'^ /It?' = eip^- It is easy to see that eq. ( ^l[ ) at most only can have solutions for negative e^- Instead of eq. (|52 
we then at large p obtain 



Va 



iV-(wo + tx)/^x , Kao « Rq\/{vq + ex)/2 



Expanding eq. ( pT| ) to lowest order in 1/p we find (instead of eq. (|53|)) 



tan 



(56) 



(57) 



which is the eigenvalue equation for the two-body subsystem with an energy corresponding to e^. If two-body bound 
states exist Cx must be one of the energies. In other words for each bound two-body state one A— value diverges 
paraboHcally according to eq. (p6|). 

For general short-range potentials, which vanish at distances larger than a given finite radius Rq, we obtain again 
the solutions described in this subsection. The Efimov conditions in particular therefore remain the same for such 
potentials. 



E. Large-distance behavior of P„„' and Q„ 



According to eqs. (||) and (|^) we must find the derivative of the Schrodinger wave function, $„ = ^ ' with 

respect to p. In the present symmetric case the three Faddeev components are equal. We choose to work in the 
first Jacobi coordinate system and must then express the orther two components in this set of coordinates. The 
transformation is given in eq. (|0) and the angular wave function is 



$ cx (f)niai) + 2i?i2 



sin(2Q2) 



4 /•7r/2-|7r/6-Qi| 

kiiai) + —7= (t>n{p,a.2)da2 , 

VO J|Tr/3-Qi| 



(58) 



where the factor 2 arises since the i?i3-operation gives the same result as that of R12. The wave function 0„(ai) 
explicitly is shown in eqs. (^ and ( ^ ) with obtained from eq. ( ^9| ) as 



iV( V A, K, ao) = 



Ah 



1 



sin(Q;oK) 



= sin('\/X 7r/6) sin(ao'\/X) — sin f (ao — it /2)^/\ 



(59) 



It is then clear that only depends on p through v A„, k„ and ao. Thus 



A$ _ K d<^n ^ dun d<^n ^ dap d<Pn 
dp " dp d\/\i dp dun dp dao 



For Pnn' we immediately obtain that 



Pnn' ^ 



duQ 

~d^ 



r'^ d 

Jo dao 



(60) 



(61) 



as seen from the orthogonality of the different Furthermore, the derivative of the overall normalization factor 
All, which also depends on uq, does not contribute to eq. ( |6T|) again due to the orthogonality of the wave functions. 
The large-distance behavior of da^ /dp is seeen from the definition in eq. (|30|) to be cx 1/p^ . 
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We now divide the integration interval into the four parts appropriate for the combination of the transformation 
between Jacobi coordinates contained in eq. ( ^ ) and the different wave functions given in the two intervals of eq. (|4^), 
i.e. 

h = [0, q;o] , h = [ao, Ti'/S - ao] , h = [tt/S - ao, 7r/3 + ao] , h = [tt/S + ao, 7r/2] . (62) 



As seen from eqs. (^1]) and (^8|) and fig. 3, the wave function in eq. ( |6l| ) is apart from Ajj independent of ao in the 
intervals I2 and /4, where the contribution to therefore is zero. 

The third interval vanishes with cx 1/ p and the ao-dependence of the wave function (again apart from Ajj) arises 
entirely from the Ah-tcrm of the two transformed Faddeev components. The derivative is proportional to Ko^dN / 9ao, 
which vanishes as ao oc 1/ p. The contribution to from the third interval is then vanishing at least as fast as 
1//. 

The size of the first interval Ii vanishes as ao The ao-dependence of the wave function arises (apart from 

All) entirely from the A/j-term in the first Faddeev component. The derivative of N is proportional to /tao, which 
remains finite together with the total wave function. Therefore Pnn' in eq. (|6^) vanishes as \/p^ for all n, n'. The 
leading order term comes entirely from Ii , where the transformed Faddeev components combined with the derivative 
of the Ah-tam of the first component contribute to the order 

For the non-diagonal Qnn' we also in analogy with eq. (|6l|) obtain 



where both terms therefore approach zero at least as The first term due to [da^/dpY oc l/p^ and the last term 

due to d'^aojdf? (xl/ p^ and the integral found in eq. (^ vanishing as \/ p. 

The diagonal term Q„„ can be written as a sum of six terms corresonding to eq. ( |63| ) and the analogous terms 

where ao is replaced by and k„. The three terms with first order derivatives of the wave function all vanish due 
to the normalization of $„. Since (dao/dp)^ oc l/p*, we are left with the leading terms 

/ da,<^*^^.^r.+ [-X^\ / rfai$;^^<i>„, (64) 



dp J Jo 

which for n ^ n' would vanish due to orthogonality. The factor (dKn/dp)^ approaches a constant when p 00 and 
so does {d\/\'i/dp)^ when v^A^ corresponds to a bound two-body state and otherwise {d\/x^/dp)^ oc l/p**. 

We again divide into the appropriate four intervals Ii, l2, 13, Ii- Let us first study the case without two-body bound 
states where the leading order term arises from the K-derivation, since then (9-\/X^/9p)^ oc l/p'*. The second and 
fourth intervals contain no K-dependence and the contribution consequently vanish. The third interval contributes at 
the most to the order i.e. ao oc 1/p from the interval size, ao oc from the second derivative of the Ah-teim 

of the two transformed Faddeev components and finally ao from the interval size of the transformation integral. In 
the first interval we have ao oc 1/p from the size of the interval and a§ oc 1/ p^ from the second derivative of the 
Ah-term in the first Faddeev component. The leading term of the total wave function arises from the two tranformed 
Faddeev components and it approaches a finite constant for p — > 00. 

The K„-dependence of the normalization constant An did not contribute to Pnn' due to orthogonality of the wave 
functions. However, for Qnn we must consider this dependence. If we define $„ = An^n, it is easy to show that 

8Arr P^^^ fl 

-^ = -\An? — $„. (65) 

Ul^n Jo >Jl^n 

Going through the different intervals we find again vanishing contributions from I2 and /4, a contribution proportional 
to from /3 and the leading order term oc 1//?^ from Ii arising from the two transformed Faddeev components 
combined with the derivative of the first Faddeev component. We therefore have dAn/dnn oc Furthermore, 
since 

da,Kg^^^n = 2Aii^J^ da.K^^^n 
^A//— / dai-l;$„ + |A,,|2 / dai$;-^l>„ (66) 

Jo Jo 
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the terms with derivatives of Ajj approach zero at least as fast as l/p^ and so did the last term as discussed above. 
Thus Qnn approaches zero at least as fast as l/p"^. 

In the case where A„ corresponds to a bound two-body state, we get the same results for the same reasons as 
described above. However, then (d\/x^/dp)^ approaches a constant for p — s- cx) and the second term in eq. (^) must 
also be considered. Although possible this is a rather elaborate procedure due to the many terms containing A„. We 
shall instead give arguments based directly on the equation of motion. 

At large distances the Faddeev components decouple as for example seen from cq. where the potential in the 
last term vanishes unless a approaches zero at least as fast as However, when a the size of the integration 
interval vanishes together with the integral itself. Thus the angular wave function is in this limit determined by an 
ordinary Schrodinger equation, i.e. 



A 

7 



(z) = 



(67) 



where we introduced the new variable z = psina « pa. The boundary conditions <j){z =—f 0) = (j){z oo) = make 
eq. ( |67| ) equivalent to the Schrodinger equation for one of the (three identical) two-body subsystems, i.e. 



U(Z) = ^U(Z) 



(68) 



where z = ry^m/ p, r is the relative coordinate, ex is the energy, p the reduced mass of the two-body system. 
From eqs. (^7|) and ( |6^ ) we obtain 



A = A + 4 = 



dzu{z) 



du 
dz 



dz^ 



(69) 



and 



y/pu{z). Introducing 4'{z) in the (5„„-equation we get 



Qnn = ^ ^ ' 



du 
dz 



dz^ 



(70) 



The integrals in eqs. ( |69| ) and (|7^) cancel in the combination (A + 15/4)/p^ — Qnn found in the coupled set of radial 
equations, see eq. (Q). The result, 2m£x/Ti^, restores the proper two-body radial asymptotics which describes the 
motion of one of the particles against the bound system of the other two particles. 

We can conclude that the coupled set of radial equations in eq. (0) decouple in the limit of large p, since both 
coupling terms Pnn' and Qnii' approach zero at least as fast as 1/p^, i.e. faster than that of the leading centrifugal 
barrier term 1/p^. This behavior at large distances is a general result valid for all short-range potentials. 



IV. THREE NON-IDENTICAL PARTICLES FOR SPIN INDEPENDENT INTERACTIONS 

In the asymmetric case with spin independent interactions the three Faddeev equations in eq. (|l6|) can by use of 
eqs. (^, (|l|), (|l|) and (^ explicitly be written as 

- ^'^r-^ + {pMp^i^^^) - A(p))(/)W(p,a.) + p'v.ipsina,) (71) 

where i = 1,2,3 and where we again omitted the label "n" . We shall follow the same procedure as in the symmetric 
case and mostly use square-well potentials, i.e. 

V^{r) = -V^'^&ir < R,) (72) 

and the reduced potentials in eq. ( [7l| ) are therefore the other stepfunctions 

v,{x) = -v^'^Q{x < X,) , (73) 
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where Wg*"* — 2mVQ''' /h'' and Xi — Rifijk- The angular Faddeev equation in eq. ( |7l| ) are now solved analytically in 
different intervals corresponding to increasing values of p. The decisive quantity for each of the square-well potentials 
is psinofi, which determines whether Vi(x) is finite or vanishes. When p < Xi we have Vi{x) 
a, e [0, 7r/2] 
Without : 
analogy to eq. (BOl) we define 



for all values of 



Without loss of generality we can assume that Xi < X2 < X3 which implies that ag^"* < ag^"* < ag'^'' where in 



= arcsin(Xi/p) 



(74) 



We shall in this section confine ourselves to small or large distances. The intermediate distances can be solved in 
analogy to the symmetric case by division into p-intervals. The related eigenfunctions arc combinations of simple 
functions, but they are rather tedious to write down in their full length and less interesting as well. 



A. Short-distance behavior: < p < Xi 

The potentials are now all constants for all a-values and eq. (^Tj) is an inhomogeneous differential equation for 
4>^''\ where the inhomogeneous part contains 0^-'^ and 0^^^^. It can therefore be solved by adding one inhomogeneous 
solution 



0^'^ {p, at) = flj sin(2naj) 
to the complete set of homogeneous solutions 



cj)^''\p,a^) = biSm{aiK,i) , Hi \f Wg'V^ + Mp) 



(75) 



(76) 



(1) 



,(2) 



eq. ([76|) can only be solutions if Vq 
section 3. Any asymmetry therefore leads to 61 
the three equations 



,(3) 



where we already imposed the boundary condition 0^'^ (p, ai — 0) — 0. 

The inhomogeneous components in eq. (^) must all be proportional to sin(2?ia), since the integrated values in 
eq. ( |7l| ) otherwise do not return the same function as required by this equation. The homogeneous components in 

and ifi = ip2 = f3 = /'i which is the symmetric case discussed in 
b2 = — 0. A set of solutions is therefore obtained, if and only if 

(77) 



p'^i'^Q^ - e)ai = p^v\^\akdj + a^dk) 
are fulfilled when k} are the three permutations of {1, 2, 3} and where we have defined 

_ sin(2n(^j) 



4n^ ~ p^e 



n sin(2(pi) 



(78) 



In the limiting case of p = we immediately find the three times degenerate solution A = An'^. 
Non-trivial solutions for p ^ only occur when the corresponding determinant vanishes, i.e. 



which determines the possible A(p). 
In the symmetric case (wg^"* = Wg^^ 



- „(3) - 



_ Jl) 


vi'^d. 


vi'U2 




'0 "3 






= 


'0 "2 


v^'^d. 


(3) 





(79) 



g = vo, di — d2 ^ d^ — do) we obtain the two solutions 
A(p) ^An^ - p^vo{l + do) 



(80) 



~X{p)^4n^-p\o{l-2do) 



(81) 



where the first is two times degenerate and the last corresponds to the symmetric solution as seen from eq. (77). 
In the general case we define 
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So = v^^'>v^^'>vl^\2did2d3 + dl + dl + dl-l) 
S, ^ vl,'\l^\l - dl) + v'i\f\l - dl) + v^^\f\l - dl) (82) 

and rewrite eq. ( [79| ) as 

+ S-ae^ + Sie + SQ = , (83) 

which is independent of p. Thus the short-distance behavior of A(p) is quadratic in p and given by eq. ([78|). 

For n = 1 we have di = —1, Sq = Si = Q and consequently the solutions e — —S2 and the doubly degenerate 
spurious solution e = 0. This is also seen from eq. ( |80| ) with n — 1. The corresponding wave function is determined 
by eq. (75) with c>c v^'^ as seen from eq. (|7 



For n = 2 we have 



, , tan^ <v5i — 1 ,„ 
cos(2^,)^ /\ (84) 
tan ipi + 1 



and 5*0 = 0. Consequently we find the spurious solution e = and the two solutions 



1 



e = - [-S2 ± ^Sl - 4^1 j , (85) 

where e is real, since 5"! > 45i for all masses and all values of the strength parameters vf' . The related wave functions 
are found from eq. ([tT]). 

For arbitrary n- values exceeding 2 we have three real solutions to eq. (^) as discussed in appendix D. For n ^ 00, 
we have — > 0. 

Since the potentials also as in the symmetric case only enter in combination with p^ , a perturbative solution to first 
order in p^ is obtained for arbitrary potentials by using Vi{x) w Ui(0) = — Wq*^- Thus the solutions in this subsection 
are also solutions for any potential in first order perturbation theory provided the depth of the square-well potential 
is replaced by the value Vi{Q) of the potential at p = 0. 



B. Large-distance behavior 

We define in this connection large distances to mean that 

(1) ^ (2) ^ (3) ^ I (1)1 ^ I (2)1 ^ I (3)1 /Q«N 

for all sets 1,2,3, which implies that a^^ < ^Mm{ipi) for all sets 1,2,3. The potentials all vanish for > and 
the solutions to eq. ( [Zl] ) are therefore the same as in eq. (^l|), i.e. 



(a,) = a, sin (^(a, - 7r/2)\/Tj . (87) 



(i) 

In the other region where Ui < we also obtain decoupled solutions, i.e. 



^'^*^(a,) = biSm{aiKi) + sin (^a^y^ 



where the first term is the solution to the homogeneous equation. The constants 6,; are therefore completely free, since 
the functions to be integrated in eq. (|7|) only involve the functions in eq. (|8^), i.e. values of > ag*'', i=l,2,3. 
The connection between Ui and Ci are now found by insertion of eqs. (|87| ) and (|8|) into eq. (^. The result is 

_ _ 2F3 

where 



= iaj+dk), (89) 
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sin ((yj, -77/2)^1) 



sin(2(^i) 



(90) 



This determines from given values of a^. 

The matching conditions at ciQ \ o/'q^ and o'q^ provide 6 Unear constraints between the parameters and bi, where 
Ci is found from eqs. (|8^)-(|90|). They only have non-trivial solutions when the corresponding determinant vanishes. 
In appendix E is shown that this condition, which determines the eigenvalues A, can be formulated as 



(91) 



where Ai and -B, are defined in appendix E. 

In the symmetric case, {Bi — B2 = B3 — Bq, Ai — A2 = A3 ^ Ao, /i = /2 = /s), we have D = {Aq — Bq)'^{2Ao 
Bq), where the solutions to = correspond to those of eqs. ( pO[ ) and (|8l|). 

Eq. ( |9l| ) simplifies in the limit of very large distances, where 



P 



(92) 



where we assumed that A remains finite. Expansion to lowest order in l/p then gives, see eqs. (E^ and (eQ), 

A, « -2F^i,, \f^a'^:L cos{X, J^) , (93) 



-^cos(X.V^) 



Bi « cos(X^\/ v^'')\/ Vq^ (^psin(^yi) + a'•^}^^^ljky^cos{^\/\)j 
where the scattering length for the two-body system in analogy with eq. (p4|) is given by 



(94) 



1 



P-jk^ 



tan X 



(95) 



The eigenvalues for p — > 00 therefore approach solutions to sin(^vA) = 0, i.e. the hyperspherical spectrum of 
A = 4n^. The equation to the next order in l/p is instead 



,7r 

tani— V A) 
^2 ' 



p 



E{i) 
'^scatPjk 



(96) 



=1 

which immediately generalizes the symmetric case in eq. ( |55| ) into 

4 ^ 

A«4n2(l-— y^ai^A^,,). 



7rp 



(97) 



i=l 



The different patological cases of extremely large scattering lengths are very different. Without loss of generality 

and consider various cases. Using eqs. (|9|), (^Sh and (|94|) we obtain 



we can here assume that I a 
of/ 

,(2) 



(1) 



< la 



(2) 



< la 



(3) 



for large values of p that the angular eigenvalues when \a'scat \ *C p <C \a!'scat \ still are given by sin(^-\/A) = and when 



kLatl < P < \ascat \ instead are given by 

± ■\/Xcos(^'\/X) sin(2iy9i) — 2sin {^Lpi — tt/2 

and finally when p <C lament | we find 

/ \/Acos(f \/a)\^ 



(98) 



VAcos(§^/A)^ (f^ + fi + fl) 



2F 



2F 



F2 



(99) 
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The last equation has the symmetric solution, \/Acos(-|'\/a) — ~4F, when all the quantities are independent of i. 
These results reflect the fact that the Efimov anomaly only is present in the last two cases where at least two of the 
scattering lengths are infinitly large with the resulting infinitly many bound three-body states obtained from eq. (|^). 
Clearly we are able to extend these expressions to any order of correction in 1/ p. 

When A/p^ diverges we can as in the symmetric case again define in the large-distance limit by A — > 2mexP^ /h^ = 
CxP^ ■ Eq. can at most have solutions for negative e^. At large p we have approximately 




ex)/ei , ^^af « X,\^ vf + . (100) 



The dominating term in the determinat in eq. ( pi| ) is B\B2B^ which diverges as p^ . To leading order in p we therefore 
obtain the eigenvalue solutions from Bi — 0. This is equivalent (see eqs. (e||) and (e|^)) to da = which immediately 
to leading order is seen to be identical to eq. (^ for the i'th two-body subsystem. Thus we obtain the general result 
that there is a one to one correspondence between the total number of two-body bound states in all the subsystems 
and the parabolically diverging A— values at large p. 

The behavior of Pnn' and Qnn' at large distance is qualitatively the same as for three identical particles. This 
result, as well as the various solutions described in this subsection, is valid in general for all short-range potentials. 
In particular, the Efimov conditions remain the same for such potentials. 



V. TWO IDENTICAL SPIN-1/2 PARTICLES 

We consider a system of two identical spin-1/2 particles (labeled "f" and refered to as fcrmions) and a third particle 
(labeled by "c" and called the core) of spin Sc- The total angular momentum J can then take the values J = Sc, Sc± 1, 
since we only consider vanishing orbital angular momentum. The spin wave functions xi*^ from eq. ( p^ ) are Xq\ 
X's ±1/2' ^"^^ ±1/2' '^t'^ich arc defined as in subsection 2.2. We have labeled the core as particle 1 and the two 

fermions as particles 2 and 3. The spin state x^i^ is not possible due to the required antisymmetrization of the wave 
function under exchange of the fermions. Since we only include s- waves the two fermions can not be coupled to spin 
1. The total angular momentum is therefore confined to be J = Sc- 

The form of the Faddeev equations depends on the spin dependence of the interactions. The fermion-fermion 
interaction, Vff{r)^ is only effective in relative s-states and no spin dependence is needed. For the fermion-core 
interaction we take the spin-dependence to be of the form 

V2{r) = V^{r) = (1 + 7,s, • Sf)Vf,{r) . (101) 

This potential is diagonal in spin space, and the diagonal matrix elements for « = 2, 3 are: 

(xi:Vi/2l^^Mlxi:Vi/2> = W = (1 + 7sy )^/c(r) (102) 
(xi:li/,|V^.(r)lxi:li/2> ^ V7-(r) = (1 - 2lif|±ll)F,,(.) . (103) 
The spin-overlaps defined in eq. (|l^) are now found to be 



/^12 


-1/2 


~ '-'0,Se-l/2 — 




+ 1/2 


_ ^13 _ 
~ "^0,^^ + 1/2 ~ 




l/2,s^-l/2 


_ /-<23 _ 
~ ^ s^ + l/2,s^ + l/2 


^23_ 


l/2,s<: + l/2 


_ /-<23 _ 
^Sc + l/2,Sc-l/2 



2Se + l ' 



— / ^c + l 

2sc + l ' 
1 

2Sa + l ' 
^4Se(^c + l) 
2Sc + l 



(104) 



The Pauli principle requires that the solutions are antisymmetric in a simultaneous interchange of all the coordinates 
of the two fermions labeled 2 and 3. This means that exchange of a2 and and exchange of the order of the couplings 
in the spin functions must give a change of sign of the total wave function. Imposing this constraint the components 
of the three-body wave function in eq. (Rlj) must be related as 
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x(3) _ - (2) 



x(3) 



.(2) 



(105) 



The corresponding Faddeev equations obtained from eq. ( |16D after integrating away the spin degrees of freedom are 
then given by 



da'i 



CHp) - p^w//(psinai)) 0^ '{p,ai) = -p^w//(psinai) 



7r/2-|7r/2-ip-Qi| , ,(2) . 

^«20;^Vl/2(P,«2, 



I /-<12 r'r/z-|i 
"•"^0,^0 + 1/2 J|ip-Qi| 

1 / ^12 ^7r/2-|7r/2-0-Qi| 

m 

'-'0,Se + l/2 J|0-Qi| 



(106) 



c'aa^f Li/2(P: "3) 



7r/2-|7r/2-i^-ai| , ,(3) / x 

rfa3?!'^/+i/2(P,a3) 



~ (A(/3) -P^i'jc(psina2)) 0lf:pi/2(P,"2) 

^«J,(psina2) [c-^^,/2^i;;/rjf/^----l rfai0^^^(p,«i) + 

(^^i%i/2,.-i/2 /r/rjr'"'""' '^«3<^i!L,/2(p, «3) 



(107) 



1 

sin(2c^) 



+ <-'s,Tl/2,Sc + l/2 JlvS-a^l «"30,^ + i/2(,P, "s) 

where the reduced potentials, w// = 2mVff/h^,vJ^ — ^raV^Jf^ ^ 
equations with second derivative Qnalrvo-rvnc! +n tVinop of en 

the constraints in eq. (105). We 



1 (JVXU\^(J^-l I^IJU^IILICHO, t^JJ ^ J J I ^ fc ^ ' ' J Cf ' ' ' 

with second derivative of a^, analogous to those of eq. (107), turn out to 
aints in eq. (105). We have then three independent Faddeev equatio 
The angles ip and ip obtained from eqs. (12) and ( |l3| ) are given explicitly by 



^, are defined as in the previous sections. The two 
/hn^\ i- i-„ identical to eq. (107) due to 



if — arctan 



M + 2m 
M 



1/2 



(108) 



= arctan 



/Af(M + 2TO)y/^ 



(109) 



where M — nii is the mass of the core and m = m2 — the mass of each of the two fermions. 

We shall again use the square-well potentials for the radial shapes of the two-body potentials between both fermion 
and core and between the two fermions, i.e. 



Vffir) - ~V,^"^e{r < Rff) , 
Vfcir) = ~V,^^'^e{r < Rf,) . 

The related reduced potentials are then defined as 

and the depth parameters for the two fermion-core relative spin states are 



-% (1 + 7.^) 



(110) 
(111) 

(112) 
(113) 

(114) 
(115) 



where Xfc = ^^2^13 and Xff = i?i/i23 are defined as in section 4. 

We shall also in this section restrict ourselves to small or large distances and omit the lengthy but straightforward 
calculations and expressions for the less interesting intermediate distances. 
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A. Short-distance behavior: < p < Min{X ff,X 



The potentials are again constants for all a-values and, as in subsection 4.1, cqs. ( |l06| ) and (107) have the general 
solutions 



(116) 



where the constants ao,i, asc±i/2,2j and asc±i/2,3 are found by substitution of cq. (HE) into eqs. (106) and ( |107| ). We 
obtain the following system of linear equations for 



(117) 



„(//) 



'^^0,Se-l/2("sc-l/2,2 + ^^^-1/2,3) + <^C'o,Se + l/2('^So + l/2,2 — 0,Sa + l/2,z) 



Se=Fl/2,2 



(118) 



^^O.SeTl/2^0.1 + ^(^sfTl/2,Sc-l/2"sc-l/2,3 + C'sfTl/2, 5^ + 1/2^*0 + 1/2, 3 ) 



where we in analogy with eq. (|7^) define 

A(p) = 



1 sm(2nip) 



1 sin(2n(^) 



n sin(2iy9) ' n sin(2(p) 

The antisymmetry expressed in eq. ( |105[) relate the unkown constants by 



asc-i/2,3 



-1/2,2 



c + 1/2,3 — ^aSe + 1/2,2 , 



(119) 



(120) 



which by use in eqs. (117) and (118) reduces the number of equations and the number of unknown constants to three. 

The solutions for e or A(p) are now obtained by demanding non- vanishing solutions for the constants Ug^i. This 
amounts effectively to three linear equations and the corresponding determinant must vanish. Thus e is obtained from 



(//) 



"'-'0,Sc-l/2 ^ 
,,(/c+) .^12 
^0 "'-'0,Sc + l/2 



^^0 ""-"0,3,-1/2 

U "'-'s,-l/2,Sc-l/2'' 



-V, 



(/c+)j^23 







Se-l/2,Sc + l/2 



ZVq ""-"o,s<: + 1/2 
Vq "'-"s,-1/2.s, + 1/2 
4^^+)(l + <f_,/,,,^„,/,) 



(121) 



This means that there are at most three A-solutions for each n; sometimes less than three, since some of them can be 
the trivial spurious solutions characterized by e = 0. 

When the determinant vanishes for e = 0, we find that these spurious solutions must satisfy one or both of the 
following two conditions: 



sin(4n93) 
sin(4(/5) 







(122) 



9 sin(4n(i?) 
n + .\. n - 



sin(4(y9) 



^sin^(2niy9) 
' sin2(2v5) 



. 



(123) 



For ri = 1 [d = d = —1) these two equations are both satisfied, meaning that two of the solutions are spurious and 
only one antisymmetric solution exists. Solving the determinant (121) the antisymmetric solution is found to be 



\{p) - \{p) + 4 = 4 - p2 



,(//) 



2Sr 



Afc+) 



2Sr 



2Sr_ 



2Sr 



(124) 



which reduces to eq. ( y7[ ) for identical and spin independent potentials 
For n = 2, where 2a^ = 1 — d, only the condition in eq. (123) is s 



satisfied, and two non-spurious antisymmetric 



solutions appear. They are the solutions of the second order equation 
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2Sr 



(125) 



+ (1 + d) 



2s, + 1+"^° 2s, 



- V, 







"0 



(1 - d) 



which combined with cq.(119) results in two solutions for A(p). 

For n > 3 none of the conditions ( |122| ) and ( |123| ) are satisfied, and three non-spurious antisymmetric solutions are 
found. For very big values of n both d and d approach zero, and the three solutions of the determinant converges 



iff) 

towards e — Vq \ e 



and e — V, 



(fc+) 



As in eqs.(pO|) and (|74|), we define 



B. Large-distance behavior 



ai"'' = &TCsm{Xff/p) , (126) 
a'o^"^ = arcsin(X/c/p) (127) 



such that the potential is non- vanishing only when the corresponding a < ao- 

Since q;q and approach zero for increasing p, we can define large distances by 

Aff) / Afc) ^ , ^ (ff), ^ , ^ Afc) 



ai"^<alf^U\^-alff^\<\0-alf^^\, (129) 



where we assumed that X-ff < Xfc- 



nil 




nil 



FIG. 4. The regions, defined in subsection 5.2, in (q/c, a//)— space arising for large distances for two identical spin-1/2 
particles. The shaded area does not enter in the present computation. 
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In fig. 4 we plot a/c = ^2 versus aff — ai where only the region, aff < Ufc, is relevant for our calculation. In the 
plot we can distinguish four different zones: 

// : a/c > Uq^^ and a// > a\^^^ V7c(/5sina/c) = and V//(psina//) = 

B : a/c > a\l^^ and a// < a\^^^ V/c(psina/c) — and V//(psina//) ^ 

C : a/c < ckg''^'^^ and a// > a\l^^ y/c(psina/c) 7^ and V//(psina//) = 

D : a/c < ckg''^'^^ and a// < ap"'^''^^ Vfdpsinafc) ^ and ^//(psina//) ^ 

Since a^ '^^ and a'^J'^'' are close to zero, regions B, C, and I? are very small. Due to the inequalities in eqs.( |l28| ) 
and ( |129| ), the integrations appearing in the Faddeev equations involve only the functions in the region where the 
potentials vanish. 

The wave functions must vanish at a = and a — tt/2 and the solutions in region II, where all potentials are zero, 
are therefore of the form 



0l^)(a,) = 6.,»sin((a, -^/2)Va) . (130) 

In the other regions (B,C and D) the same form of the wave function is a solution when the corresponding potentials 
vanish, i.e. for a^ > ag*^ . When a^ < ag ^ in these regions, we have solutions of the form 

0^'^(ai) = Cs,i sin(aiKs) + ds,i sin (^UiV^ , (131) 

where the first term is the solution to the homogeneous equation and consequently the different Kg-values are given 
by 



The constants Cs_,; are therefore completely free, since the functions to be integrated in eqs.( |l06|) and ( |107| ) only involve 
the functions in eq.(130), i.e. values of ai > ag*' for all three components, see fig. 4. The constants bs.i and ds^i are 



linearly related through the Faddeev equations. 

The detailed solutions in the different regions are given in appendix F. The eigenvalue equation again takes the 
form of a vanishing determinant where the matrix elements are given in eq. (F^3|) . 

In the limit of very large distances we can for finite A— values use an expansion to lowest order in 1/p as given in 
appendix F, where also the case of diverging A— values are considered. The linear set of equations in eq.(F|l2|) then 
reduces to 



A, fyicos(;yi)^/i23 +sin(jyX)^ = -aoVl^ 

y 2 p ^ J P 



iff) 

'^A^23 (133) 
P 



A2 V^cos( Jv^)^^A^i2 + sin(Jv^) = -a,^_^/^Vl^^^ii^2 (134) 
\ 2 p 2 / ' p 



( (/c+) \ (/c+) 

A3 yicos(5v^)^^/.i2 + sin( Jv^) = -a,^+i/2V^^^/ 

\ 2 p 2 y ' p 



7^12 , (135) 



where /X12 and p^z are constants defined in appendix A, and oq, as^±i/2 depend linearly on the coefficients Ai as defined 
in appendix F. The eigenvalues for p — 00 therefore for finite scattering lengths approach solutions to sm{J^\f\) — 0, 
i.e. the hyperspherical spectrum of A = At? . 

The solutions to the next order in 1/p only exist when the corresponding determinant is zero which then defines 
the behavior of A(p) for large values of p. We have previously assumed that ag'^''^^ < ag'^"-' and both quantities are 
close to zero in the large-distance limit. If we further assume that ag — or equivalently Xff — 0, only regions II 
and C survive in fig. 4. Then Ai = 0, \ascat\/P ~ i^'^^ finite aj^^l and large p) and the fermion-fermion interaction 
completely disappears. The determinant is reduced to a 2 by 2 determinant, that gives the following expression for 

Hp) 
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Afc+) (fc-) 
scat ^scat 



2sc + l 

-A cos^ 



\0'scat 



Afc-) 



Afc-) 



(136) 



M?2 - 4Ai!2/' = 



where / and / are defined in eq. (fH) . 

In the Umit where both scattering lengths are large compared to p, i.e. |ai{af^| S> p, we arrive again at eq. 
The occurrence conditions and properties of the Efimov states can be derived from eq. (|l36| ). The details are discussed 
in Hi. 

The behavior of P„„' and Qnn' at large distance is qualitatively the same as for three spinless particles. This result 
is again, along with the various solutions described in this subsection, valid in general for all short-range potentials. 
In particular, also the Efimov conditions remain the same for such potentials. 



VI. NUMERICAL ILLUSTRATIONS AND GENERALIZATIONS 

The analytical solutions in the previous sections are derived for square-well potentials. The results are similar for 
other short-range potentials. Illustrations by use of smoother gaussian potential are therefore appropriate and we 
shall first show numerically calculated angular eigenvalue spectra for different symmetries of the three-body system. 
Then we shall give a survey of the analytical results for different distances and show how the results can be used 
for arbitrary short-range potentials. The simplest case of intrinsic spins of the particles were also considered in the 
previous sections. Other cases with more complicated spin structures can be worked out. It is straightforward, but 
results easily in rather extended formulae. We shall here indicate how to proceed in the case of two different spin-1/2 
particles. 
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FIG. 5. The lowest angular eigenvalues A = A — 4, when all orbital angular momenta are zero, are shown as function of p/b 
for the case of three identical bosons for a gaussian potential, Voh^ /mh^ exp(-(r/6)^), where m is the mass of the bosons. The 
actual value of Vo = —3.08 corresponds to a bound two-body state. 
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A. Numerical illustrations 



The angular eigenvalues, still only s-states, for the symmetric case of three identical bosons are shown in fig. 5 
as function of p/h for a gaussian potential. The parameters of the potential correspond to a two-body bound state. 
The lowest eigenvalue consequently bends over and diverges parabolically as described in eg. (p6|) . The higher lying 
levels then come down and the hyperspherical spectrum seen at p = is approached at large distances. The values of 
A = in? — 4 at p = are 0, 32, 60 and the spurious state starting from A = 12 corresponding to n = 2 in eg. (p9|) is 
omitted. Due to the symmetry reguirement only one state appears for each value of n. 
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FIG. 6. The lowest angular eigenvalues A = A — 4, when all orbital angular momenta are zero, are shown as function of p/b for 



three different spinless particles for gaussian potentials. 



V^^n'^/mih' exp(-(r/6) = 



The actual values Vr 



-1.28,-1.54, 



and rrii = mi , 2mi , 3mi , respectively for i 
no bound states in the other subsystems. 



1, 2, 3, correspond to a bound two-body state in the subsystem labeled 



-3.08 
1 and 



In fig. 6 we show the angular eigenvalues as function of p/b, again only s-states, for the asymmetric case of three 
different spinless particles interacting via gaussian potentials. Only one of the potentials has a bound two-body state 
and conseguently the lowest eigenvalue bends over and diverges parabolically as described in eg. (100). Again the 
higher lying levels come down and the spectrum at p = emerges also in this case at large distances. The spurious 
level corresponding to n=2 is still omitted, but now two other levels appear at A = 12 corresponding to (spatially) 
non-symmetric configurations. Also for n=3 at A = 32 two more levels corresponding to non-symmetric states appear 
in addition to the totally symmetric state. For all values of n > 3 the structure of one symmetric and two asymmetric 
states remains unchanged. 

In fig. 7 we show analogous numerical results for two identical spin-1/2 fermions plus one third particle (core) 
of spin Sc = 3/2. Now one additional degree of freedom appears. It is related to the two different fermion-core 
spin couplings. However, simultaneously the two fermions are restricted to totally antisymmetric two-body states. 
This clearly removes some of the otherwise possible states. Since the hyperspherical spectrum always both is the 
starting point at p = and the asymptotic limit for large p, the combined result is a spectrum very similar to that 
of fig. 6 corresponding to three different particles without spin degrees of freedom. If the core had been spinless, the 
fermion-core spin coupling had been unigue and one state less would have appeared for every value of n > 2. 
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FIG. 7. The lowest angular eigenvalues A = A — 4, when all orbital angular momenta are zero, are shown as function of p/h 
for gaussian potentials, Vq^^ h"^ / mib^ exp(-(r/ti)^), for a system of two identical spin-1/2 fermions plus a third particle of spin 
Sc = 3/2. The actual values fJ*' = -6.16, -0.99, -1.09 and rm = 2m2, 1712, m2, respectively for i = 1,2,3, correspond to the 
fermion-fermion interaction, the fermion-core interaction for relative spin of 2 and the fermion-core interaction for relative spin 
of 1, respectively for i = 1, 2 and 3. 

B. Survey and generalization of the solutions 

Various analytical solutions apply to different regions in the two-dimensional (p, a) coordinate space. A survey is 
shown in fig. 8 where the regions labeled I,A,B,C and II refer to the division in the subsections of the symmetric 
case. (The more general asymmetric case requires more divisions at intermediate distances. However, the picture 
basically remains unchanged.) The simplest cases of small and large distances (I and II) actually constitute a very 
large fraction of the total coordinate space. These solutions are also appropriate for general short-range potentials 
where the approximate validity results from perturbative treatments. 

At large distances, where ao and k can be expanded to first order in we obtain considerable further simplifica- 
tion. The eigenvalue equation can in this limit be expressed entirely in terms of the scattering length of the potential. 
In fig. 9 we show the lowest angular eigenvalue for a square-well potential as function of p. We compare with the 
approximate solution at very large distances given in eq.(|5^). The 1/p asymptotic convergence towards the limit is 
clearly seen. This approximation is normally too inaccurate for the interesting shorter distances. 

A numerical procedure for general short-range potentials now suggest itself. First we construct the equivalent 
square-well potential with the same scattering length and effective range as the potential in question. The solutions 
are then obtained analytically in a— space as described in this paper. They can be anticipated to be very accurate 
solutions except in smaller regions at intermediate distances, where direct numerical integration interpolating between 
the available analytical solutions then must be used. The gain in accuracy and speed is substantial. 

An example is shown in fig. 9 where we used an attractive gaussian potential of range b and strength —2.19h^ /mb^ 
with a scattering length, ascat = 56, and an effective range, i?e// = 1.646, where b is the range of the potential. The 
corresponding radius and depth of the square-well potential is then Rq = 1.476 and Vq = — 0.93?i^/m6^, respectively. 
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intermediate 
distance 




FIG. 8. The regions in a-space as function of p/Ro for a square-well potential of radius _Ro for three identical bosons. The 
shaded areas indicate the small and large-distance regions I and II, where the wave functions are particularly simple. 

The curves for the two potentials are indistinguishable at large distances down to p = « 26. Then the 

square-well result is lowest until p ~ b, where the curves cross and the gaussian result stays as the lowest all the 
way down to zero. The largest deviation is less than about 0.75 compared to the minimum value of about -6. This 
relatively small region, 0.3 < p/b < 2, at intermediate distances, corresponds to the surface region for the two-body 
potential. The perturbative short-distance solution for the gaussian potential is accurate up to about p ~ 0.36. The 
large-distance behavior describing the approach towards the asymptotic limit is only a reasonable approximation at 
very large distances. 

The exact square- well solution is almost quantitatively a good approximation. The two energies would be close due 
to the similarities and the fact that a larger potential at smaller distance is compensated by a smaller potential at 
larger distance. The square well tends to confine the wave function in a somewhat more narrow region around the 
minimum. 

The (smaller) differences between the two potentials occur in an essential region, where a substantial part of the 
attractive pocket of the effective three-body potential is contained. High accuracy for an arbitrary potential therefore 
requires a treatment better than that corresponding to the square- well solution. In practise the most efficient procedure 
is numerical integration starting with the small-distance perturbative solution. Reaching large distances corresponding 
to p < RoV^ the exact square-well solution, which then is particularly simple, can be used to a very high accuracy. 



C. Two different spin-1/2 particles 

We consider a system consisting of two non-identical spin- 1/2 particles and a third particle of arbitrary spin Sc- 
The procedure is then the same as in section 5, but now the spin state Xi^'' is allowed and the antisymmetrization 



constraint eq.(I05) is not required. The three Faddeev equations can then be written as a system of six equations, 
where the angles ipi^ ip2, and 1^93 all are different. Since the two fermions can couple to spin 1, the total spin of the 
three-body system can be J = Sc, Sc ± 1- In case when J = Sc ± 1, the spin states Xo \ ^^'^ xf\i/2 ^'^^ 

possible and the total wave function is symmetric under exchange of the two fermions, and the number of Faddeev 
equations is again reduced to three. 
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FIG. 9. The lowest angular eigenvalue A = A — 4 as function of pjh for a system of three identical bosons. The interactions 
correspond to a gaussian potential (solid curve), V^h^ jm}? exp(-(r/&)'^), Vo — —2.19, and a square-well potential (short-dashed 
curve), radius i?o = lA7b and depth — 0.937i^/m6^, with the same scattering length as, 



The perturbative solution for the gaussian potential at small distance (eq. 
(eq. (pq)) are also shown (long-dashed curves). 



= 56, and effective range Reff ~ 1.64fe. 
J)) and the limiting behavior at large distance 



At short-distances the general solutions are given by eq.(116). They give a 6 by 6 determinant when J ~ Sc and a 
3 by 3 determinant when J = Sc ± 1- Making these determinants equal to zero and solving the equation for e = one 
can extract the spurious solutions and the non-spurious ones for e ^ 0. 

At large distances the procedure is also analogous to the one shown in the subsection 5.2. The only difference is 

„(/l/2) Afic) , (he) 



that now we have ag^''^', ckq^''', and a^^''', and the number of regions we have to consider (similar to regions II, 



A, B, and C in subsection 5.2) is much larger. Technically the problem is more cumbersome, but conceptually it is 
identical. 



VII. SUMMARY AND CONCLUSIONS 



The new method to solve the Faddeev equations in coordinate space has several advantages especially in the long- 
distance description. The reason is the access to the asymptotic analytic properties made possible by the two-step 
sequential procedure. The generalized angular equations are first solved for each average distance (radial coordinate) 
between the three particles. The periodic behavior or the finite intervals characteristic for angular variables provide 
discrete quantized solutions which subsequently are used as a complete basis set in an expansion of the total wave 
function. The expansion coefficients (as well as the basis) are functions of the radial coordinate. They are in the next 
step determined together with the total energy from the coupled set of radial equations. We have first formulated this 
general procedure for s-states with the inclusion of intrinsic spins for each of the three particles. For spin-independent 
interactions we recover the previous general s-state Faddeev equations. 

We study the symmetric case described by three equal Faddeev components. We assume spin independent interac- 
tions and arrive at one integro-differential angular Faddeev equation for s-states. This equation is solved analytically 
for a square-well potential. The angular wave functions are especially simple at small (one sine function) and large 
distances (one or two sine functions). These solutions are, with an appropriate interpretation, more general than their 
origin as square-well solutions seems to suggest. At small distance, they are the first order perturbative solutions for 
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an arbitrary potential with a value at the center equal to the square-well depth. At sufficiently large distances where 
the general short-range potential has vanished the solution obviously is identical to the square-well solutions outside 
its radius. We show explicitly that the coupled set of radial equations decouple at large distances. 

At intermediate distances the solutions are linear combinations of up to four sine functions of varying arguments. 
The corresponding eigenvalues are explicitly given at small distance where they vary parabolically with distance 
starting from the hyperspherical spectrum. At intermediate distances the eigenvalues are solutions to cubic equations 
and a simple trancendental equation involving trigonometric functions arises at large distances. 

The length scale is here defined by the radius of the square-well potential and the large-distance solution is valid 
when the average radius exceeds Rq\/2. Expansion of the eigenvalue equation to lowest order in the inverse distance 
brings additional simplification, but introduces also the two-body scattering length as another length parameter. The 
Efimov condition leading to infinitly many bound three-body states is then obtained and seen to be independent of 
the short-range potential. 

The asymmetric case, still s-states and spin independent interactions, is described by three Faddeev equations. We 
can again solve the angular part analytically for squarc-woU potentials ac;ting betweeen each pair of particles. Each 
Faddeev component of the wave functions is at small and large distances of the same form as in the symmetric case. 
These solutions are again appropriately interpreted solutions to general potentials either in perturbation theory or 
asymptotically at large distances. The corresponding eigenvalues are in this case obtained either as solutions to a cubic 
equation or to a trancendental equation involving trigonometric functions. At intermediate distances the solutions 
can still be found analytically to be combinations of exponentials and trigonometric functions. However, the number 
of terms is now substantial and wc abstain from writing down these soulutions. 

Expanding the eigenvalue equation to first order in the inverse distance again brings additional simplifications and 
introduces the three (diff'erent) scattering lengths. The Efimov conditions leading to infinitly many bound three-body 
states can then be seen to exist when at least two scattering lengths are infinitly large. Asymptotic large-distance 
behavior of the effective radial potentials are derived to first order in the inverse distance. 

The case of two identical spin- 1/2 fermions and one additional core particle of arbitrary spin can of course also be 
described by three Faddeev equations. We still include only s-states and assume a spin-spin type of interaction between 
fermion and core particles. Elimination of the spin degrees of freedom results in three angular Faddeev equations 
where the reduction from six to three equations arise due to the requirement of an antisymmetric wave function. 
Again these equations are solved for radial square- well potentials both for small and large distances. Expansion to 
first order in the inverse distance expresses the angular eigenvalue equation in terms of three (different) scattering 
lengths. The Efimov conditions are described and discussed. 

We finally considered qualitatively the case of two different spin- 1/2 particles interacting mutually and with a third 
particle. The original three Faddeev equations results now, still only including s-states, in sixcoupled angular equations 
after elimination of the spin degrees of freedom. Small and large-distance behavior can be studied analogously for 
radial square- well potentials. The systems considered here containing spin-1/2 particles are chosen as the simplest 
examples of three interacting particles with an intrinsic structure showing up in the form of spin degrees of freedom. 
Other similar examples can be worked out in the same way. 

The angular solutions obtained with square- well potentials have properties characteristic for short-range potentials. 
To calculate the total solution we must in addition solve the coupled set of radial differential equations. This is a 
rather simple numerical problem. The coupled set of radial equations decouple at large distances. Only s-waves are 
included in the discussion, but they are usually the most interesting components in the wave functions. For bound 
states only s-waves extend far beyond the radii of the short-range potentials. They decouple from higher angular 
momentum components at large distances and often constitute the dominating part of the total wave function. 

The square-well solutions are intrinsically interesting. Furthermore, they approach at large distance the solutions 
for arbitrary short-range potentials. The exact solutions are valid down to distances close to the radius of the square- 
well potentials. The solutions for general potentials can now be found by adjusting depth and radius of a square-well 
potential to obtain the same scattering length and effective range. The large-distance solutions obtained as described 
in this paper are then correct down to times the radius of the square- well potential, i.e. at much smaller distances 
than given by the large-distance expansion for general potentials. Therefore these results can be used to make the 
numerical computations both substantially faster and more acciirate. 

In conclusion, we have discussed a method and derived solutions providing basic insight into the Faddeev equations. 
Practical numerical improvements are suggested for arbitrary short-range potentials. 
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Appendix A: Coordinates 

We consider a system of three particles with masses rrii and coordinates r^. The Jacobi coordinates are defined as 



X, 



i — f^jkYjk,yi — fJ-{jk)i^{jk)i 

1 \ 1/2 / 1 ^ I ^ \ 1/2 

1 nijirik \ /I [mj + rnkjuii ^ 



fJ-jk = ,^J'(jk)i = (1) 

\m nij + mfe J \m mi + m2 + J 

TTljVj + TTLkYk 

rjfc = - rfe,r(jfc), = r, . 

" ' "^j + TOfc 

where is a cychc permutation of {1,2,3} and /i^ are the reduced masses of the subsystems in units of an 

arbitrary normalization m. 

The hyperspherical variables are introduced as 

p, ria;, =x7|xi|, YVy^ =yi/|y,|, a.,, (2) 

where at is in the interval [0,7r/2] 

p^=Xi^+yi^, |xi| = psina^, |yi|=pcosai. (3) 

We omit the indices where we need not emphasize the particular set of Jacobi coordinates. Note that p is independent 
of what set is used. 

The relation between three different sets of Jacobi coordinates is given by 

Xfc = Xi cos ifiik + Yt sin Lp^k , Vk = -Xi sin Lp^k + y* cos ipik (4) 

where the transformation angle (pik is given by the masses as 



LPik = arctan (5) 

where (—1)^ is the parity of the permutation {i, k,j}. 

Appendix B: B-region solution at intermediate distances 

If the wave function in eq. (p&) is a solution to eq. (^5|) we must have that 



^//g.VIV3 + -±=B'J = ^lle^Vl^/^ (1) 
iV3X V3A 



iV3A V3A 



(4) 



for each of the values of k=l,2,3 and furthermore 



^11^ ( gi%/I(Qo-7r/2) _|_ g-iV^(ao-7r/2) \ 
2\/I V / 



^'i\/I(27r/3-Qo) _ ^//g-i\/I(27r/3-ao) 



(5) 



(fe)^,j<,"'(2V3-ao) _ ^0)g-K(,"'(27r/3-ao) 



) + /2V3-ao '?^C'(p, a')da' = 0. 
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Eqs. (B|i|) and (B|) have for A 7^ 16/3 the unique solution 

1 - ' 



^//giVA7r/3 



iV3A 2Ajj 



16 
3A 



(6) 



1 



1-16 ^ 
3A 



'3A 



(7) 



For A — 16/3 the two equations determining B^^ are identical and infinitely many sets of coefficients exist. The 
constraint, in addition to eq. (B||), on are then 



4 



'3A 



Eqs. (B|) and only have non-trivial solutions when 



[k'^') +voP +A = ± 



V3 4) ' 



(8) 



(9) 



which for each of the possible signs is a cubic equation with the three complex solutions Kg where k=l,2,3. The 
corresponding coefficients are then related by 



±iB 



i^) -_Bfe-'B-/\ fc = l,2,3. 



(10) 



The solutions corresponding to the different signs are related by an interchange of the B'^jf' and i?!'''' terms in eq. ( |3' 

With the expressions in eqs. (b|), (B0) and (b|i|) for the coefiicients we can rewrite the link to the C-region from 
eq. (B|) 



as 



rao 3 r r>(*:) , 



An 1 



fc=l L '^B 



e B 



A(7r/2-ao)j-2 ^ e^'"^ A(7r/2-Qo) j-j^ 



^3A 



(11) 



Using the wave function (pc explicitly together with the expressions for the C-coefficients from appendix C, we 
calculate the integral in eq. (b||) and arrive at 



raa 3 
/ 4>c{p^a')da' 



(12) 



where the summation excluded fc = 0, since this contribution vanishes. Combined with eq. (BpT|) this gives one 
constraint between the coefficients An, B^^^ and A^^'^ for k=l,2,3. 

Appendix C: A- and C-region solutions at intermediate distances 

If the wave functions in eqs. (B^) and (|38|) are solutions to eqs. (pS) and (|3) we must have that 



V3 



(1) 



V3 



(2) 
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V3 



(3) 



for each of the values of k=0,l,2,3 and furthermore 



''/If? ^ 



(27(_'=)g-K^"^(27r/3-ao) 



(4) 



Eqs. (C|), (C|) and (C|) only have non-trivial solutions when the determinant vanishes for this linear system of 
equations. This is equivalent to either 



or instead, when eq. (CH) is false, that 



vap + X = ±i 



Avop^ 



(5) 



(6) 



The first of these equations in eq. (C|^) is a second order equation with two solutions for k[^q. They only differ by a 
sign in k and they therefore correspond to the same wave function and count as one solution here labeled by fc = 0. 
The corresponding coefficients are then related by 



(7) 



The expression in eq. (C|^) is for each of the possible signs a cubic equation with three complex solutions k[^q here 
labeled by k=l,2,3. The corresponding coefficients are given by 



t/3 



AW 



(1 ± iVs) 



(8) 



for the values of k=l,2,3. The solutions corresponding to the different signs are related by an interchange of the C\_ 
and Ci'^-' terms in eq. (^8|). 

With the expressions in eqs. (CQ) and (C||) for the coefficients we can rewrite the link to the B-region from eq, 

as 



(fc) 



Jao-V3 ^B{p,a)da 
^ ((1 ^ i/V3)e-Al-('-o-^m + (1 ± i/V3)e--A>o-m) 

An ^ ' 



(9) 



Here the summation does not include fc = 0, since this contribution vanishes. Using the wave function 0^ explicitly 
together with the expressions for the B-coefficients from appendix B, we calculate the integral in eq. (C|) and arrive 
at 



1-n IZ—OLQ 
ao — 7r/3 



— ''■/3) I iiAii 



which combined with eq. (C||) gives one constraint between the coefhcients An, B^'' and A'^''^ for k=f ,2,3. 



7r/6) 



(10) 



Appendix D: Properties of the asymmetric solutions 

We have not found a general rigorous proof for the claim that three real solutions exist for n > 3. However, 
the numerical computations all unanimously support the theorem and several very different limiting cases shall be 
discussed below. The intermediate cases are probably similar and with a determined effort possible to prove. 
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When n increases df approaches zero and 
^0- 



(1) (2) (3) a ^ ,,(1),,(2) I ,,(1),,(3) I ,/2),,(3) o 
: -Jl ^ «0 ^0 + ^0 + ^0 ' '^2 



„(1) „(2) 



(3) 



(1) 



with the solution e = Wq*' to eq. (p3[). 

In general the cubic equation in eq. (p3|) has three real solutions if 



< 



which can be rewritten as 



1 



(3) 

When one of the potentials, for example Vq , is very small or vanishes we get in this limit 

S'o , Si^ v'a 'I'o (1 - da) > S2-^ ~vy ~vy 
and the condition in eq. (D^) becomes 



1 /I 



(2)\2 , „(1)„,(2) ,2 



27 ' 108 ' 27 V4 

When all Vq^ are equal the condition in eq. (D||) reduces to 



,,(1)^2, 



27' 



27 



1 2 1 2 1 

+-y -y - -xy < 



where the overall factor (wg''')^ has been removed and the new variables x and y are defined as 

X ^ 3- dj- dj- dj , 2did2d3 + dj + dj + dj - 1 . 
Direct computation of the quantities reformulates eq. (D|^) into 



— {dl+dl+dlf + {did2ds)^<0 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



which is fulfilled for any set of di 



Appendix E: Eigenvalue equation at large distances for three different 
particles 

The matching conditions at a'^\a^Q^ and a^^^ at large distances for the asymmetric case provide the eigenvalue 



equation. This is found by equating the functions in eqs. ( |87| ) and ( |88[ ) and their first derivatives, i.e. 

ai sin ^(aQ*-* — 7r/2)\/J^ — Ci sin ^Oq''' "v/ 
ai cos ^(ofg ' — TT/2)\/Xj — a cos (^aQ^\/Xj 



= 6j sin(a|)*^K^) 
A = biKi cos^aQ^' Ki) 



(1) 
(2) 



and then eliminating bi resulting in 

((a«-7r/2)v^) 



a,- sm 



— c, sm a, 



o^Vx) cos(a[)'^Ki) - 



Gi cos 



A — Ci cos 



(3) 



Furthermore eliminating by use of eqs. (|89|)-(|90|) then provide 3 linear equations in a^. They only have non-trivial 
solutions when the corresponding determinant, D —det{dik} vanishes. The matrix elements are 

da = Ki sin (J^Uq^ — tt/2)\/x^ cos{aQ ^ ki) — \/\cos (j^ctQ^ — Tr/2)\/J^ sin(ag*''Ki) (4) 
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dik = , for i ^ k 



where fj is defined in eq. (pO|) and 



Ki sin (^aQ^\/Tj cos{aQ^ Hi) — V^cos (^aQ^\/\^ sm{aQ^ Ki) 



By further defining 



d P 



the determinant can then be written as 

D = B1B2B3 + 2A1A2A3 - B1A2A3 - A1A2B3 - A1B2A3 
The eigenvalues A are then determined by D — 0. 



(5) 



(6) 



(7) 



(8) 



Appendix F: Solutions to the case of two identical spin-1/2 particles 

The solutions to the Faddeev equations in eqs. ( 106 ) and (107) are first found independently for each of the four 
regions of fig. 4. We shall use the notation ajf — ai and ajc = 0:2- In region II , where all potentials are zero we have 

(1) 



P{a^) ^ Ai sin (^{a^-n/ 2) \/x^ , 



where (i, s) = (1, 0), (2, Sc — 1/2), (3, Sc + 1/2) and eq. (105) relate these to the remaining spin components of i/i'^^ 
and 0(2) . 

In region B of fig. 4, where only one of the potentials is identically zero, we have instead 

(2) 
(3) 



(4) 



(/'o^''(a//) = 61 sin(a//Ko) + aosin(Q;//VA 



(a/c) = 6i sin I (a/c - 7r/2) V A) , 



where (i, s) = (2, - 1/2), (3, + 1/2), = V t^^^-^'V^ + Kp) and 



4(Cd,^s, -1/2^2 - C'o.^s, +1/2-43) 



where we in analogy to eq. (b0|) define 



sm\^{ip-TT/2)VXj , sin ( ((^ - 7r/2)\/A 

f ' 771 \ ' ~ 



sm{2ip) ' sin(2(^) 

In region C of fig. 4, where the other potential is identically zero, we have 

= ci sin(|(a//-7r/2)VA 
(/)^*)(a/c) = Ci sin(a/cKs) + as sin(a/c Va) , 



where {i, s) = (2, Sc - 1/2), (3, + 1/2), Ks^±i/2 - ^v^^'^'^p^ + A(p) and 

2 

71 



^0,s, -1/2^1/ + (C's^_i/2,s^_l/2^2 + (^8^-1/2,5, + 1/2^3)/ 



flsc + 1/2 



Vx 



Cq^s, + 1/2^1/ + (~C'sf-l/2,s, + l/2^2 + (^8^+1/2,3, +1/2^3)/ 



(5) 

(6) 
(7) 

(8) 
(9) 



In region D of fig. 4, where both potentials are finite, we have 
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(10) 



where (z, s) = (1, 0), (2, - 1/2), (3, + 1/2). 

Imposing continuity of the functions and their first derivatives at the boundaries between the different regions of 
fig. 4, the constants must be related by 



d = Ai , b2 = A2 , bs ^ A3 , di = bi , d2 = C2 , ds = C3 , 
and the remaining six constants obey the following set of six linear equations: 

Ai sin ((4^^^ - = bi sinia^J^^o) + gq sin{alf^^\/x 

Ai- ^ 



[ID 



\/\cos (^{a''/-^-' — ^)\/Xj — biKQCOs{aQ^'^^ kq) + aoVJ. cos{aQ^ \/\) 



A2 sin I 



(("o''''''' - f )VA^ = C2sin{a'^J''^K,^_i/2) + a^^_i/2 sin(a[)-'''=^ \/I) 



(12) 



Ao 



VXcos (^{a^f''^ ~ f )^) = C2Ks,-i/2COs(a^-''''''Ks^_i/2) - as^_ili^c.os,{a-^''^ 
Az sin {[o^q"^ - f )V^) = C3 sm{a'^''^ K^^j^xl2) + as,+i/2 ^Ma^i''^ \f\) 



iVTcos ^(Q!o 



A = C3K, 



se+1/2 



COS 



{oli'''' i^s,^\l2) + 0^^+1/2 VAcos(a[)''^''^yA) 



Since only depend on Ai we can easily eliminate bi, 02,03 from the equations in eq. (F|l2|). This leaves 3 
homogeneous linear equations in Ai, A2, A3 with a corresponding determinant D —det{dik}. The matrix elements 
are given by 



D2 + F2zCll_i/2,s,-l/2 ' '^s 



3 - -O3 + -F32C'^f+i/2,s, + l/2 ' 

•^ii — Di , dik — FikCl'^g^ , 



where si — 0, S2 = Sc — 1/2, S3 = Sc + 1/2, the spin overlap functions are defined in eqs. ( [l5|) and 
A = Ks sin (^(a^)*^ - f''^) <=os("o'''«s) - \/Xcos (^(q;[,''' - ^)^) sin(a[)'-'Ks) , 



(13) 
(14) 

(15) 



K» sm 



cos(aQ ^Ks) — VTcos ^Ctg sin(Q;Q ^Kg) 



2(l + 5.,i) 



where fj is defined in cq. (pO|). 

The values of A are as usual determined from D = 0. 

Simplified expressions can be obtained in the limit of very large distances where 



-2(1 + (5z,i)Mjfe V w^'^ai'lt cos(X,a/4'Vj 



(16) 



(17) 



D, w - cos(x,y vy)\/vy [psm{-v A) + ay^t^fj-fe V Acos(- V a; 



(18) 



Here i = 1,2,3 corresponds to ff,fc—,fc+, respectively. Then Vq is defined in eqs. (114) and ( |115| ) and the 
scattering lengths, a[{.ll, a[{.'g^\ o!'Icat ^ fo'" the two-body systems are defined in analogy with eq. (95) in terms of 
radii, potentials and reduced masses. 

When A/jO^ remains finite the dominating terms in the determinant are Di. Thus D is approximately diagonal 
with the matrixelemnts Di given in eq. (F]!^) . The eigenvalues in the limit of large p then again corresponds to the 
two-body bound states obtained by using eq. (F^5|) and solving Di — Q. 
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